AD-A278  767 


RL-TR-93-264 
Interim  Report 
December  1993 


GREEN’S  FUNCTIONS  FOR  AN 
ANISOTROPIC  MEDIUM:  PART  I 
UNBOUNDED  CASE 


ARCON  Corporation 


Saba  Mudaliar 


DT1C 

ELECTE 
MAY  0  2  1994 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


94-12959 


DTIC  QoALili  5 


Rome  Laboratory 
Air  Force  Materiel  Command 
Griffiss  Air  Force  Base,  New  York 


94  4  28  059 


This  report  has  been  reviewed  by  the  Rome  Laboratory  Public  Affairs  Office 
(PA)  and  is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At 
NTIS  it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RL-TR-93-264  has  been  reviewed  and  is  approved  for  publication. 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  Rome  Laboratory 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  RL  (  ERAA  )  Hanscom  AFB  MA  01731.  This  will  assist  us  in  maintaining 
a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on  a 
specific  document  require  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


Puble  mpotogbkJdrtfcrMiariMlGnrfMomMlDn  to  MtnMdtoMngi  t  houi  am  wnm  tTaudng  itn  BT»  tar  Ipurg  rtrucoy  mmayrg  mamrg  amt  iojcx 

B«gwne trcnmTmrQ ty Omt n— otl  «nd og’vmj-\i gTd mn—rg tfn  nct^lr-^ d *jrnrmkx\  S«no  ccnry  t«  mErdnBBT  Ol*q»-i  — otm» or  try  arm  tno  at  ins 
colMtantf  HgWte\heftrtmi«uqgElBn»tanlLBlnB»*ai«dHaW»MnpBnHMd»—»Sf»fciOt»a»—ta  Irtan—nnOp— gn»ida«po«  izisjrft* «on 
DfllgNl.  Sute  <»■  V»  WlWOg.  1dtB»«0IWBiB>  M— pTwt  »id  Bu^t.  mrftaaueaonPrcpO  (tgOt-011,  WgdTrqmy  DC  20SCO 


1.  AGENCY  USE  ONLY  (Lmv«  Blank) 


2  REPORT  DATE 

December  1993 


4.  TITLE  AND  SUBTITLE 

GREEN'S  FUNCTIONS  FOR  AN  ANISOTROPIC  MEDIUM:  PART  I. 
UNBOUNDED  CASE 


&  AUTHOR  (S) 

Saba  Mudaliar 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS (ES) 

ARCON  Corporation 
260  Bear  Hill  Road 
Waltham  MA  02154 


a  REPORT  TYPE  AND  OATES  COVERED 
Interim  Jan  93  -  Oct  93 


a  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESSES) 

Rome  Laboratory  (ERAA) 

31  Grenier  Street 
Hanscom  AFB  MA  01731-3010 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


RL-TR-93-264 


1 1 .  SUPPLEMENTARY  NOTES 

Rome  Laboratory  Project  Engineer:  Zachary  0.  White/ERAA/(617)  377-3191 


12a  DISTRIBUTION/AVAILABIUTY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


1 12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Mtrtrvtn  200  word4 

Dyadic  Green's  Function  (DGF)  for  layered  anisotropic  media  is  essential  for  the 
electromagnetic  field  analysis  of  several  problems.  With  the  goal  of  deriving  the 
DGF  of  a  two-layer  biaxially  anisotropic  medium  we  derive  in  this  report  the  DGF  of 
a  corresponding  unbounded  problem.  Using  the  Fourier  transform  method,  an  auxiliary 
dyadic  Green's  (ADGF)  is  first  derived.-  The  DGF  is  then  obtained  by  performing  a 
simple  linear  transformation  on  the  ADGF.  It  is  expressed  in  a  compact  dyadic  form 
in  terms  of  two  characteristic  waves,  viz.,  the  a-wave  and  the  b-wave.  Some  features 
of  the  DGF  are  discussed  by  comparing  our  results  with  those  of  a  corresponding 
uniaxial  problem. 


,  DTIu  ^uALiT  X  iiiiijLi  ilCTiSD  3 


14.  SUBJECT  TERMS 

Green's  Function,  Electromagnetic  Waves,  Anisotropic  Medium 


1ft  PWCC  CODE 


1 7.  SECURITY  CLASSFICATION 
OF  REPORT 

UNCLASSIFIED 


1  a  SECURITY  CLASSFICATION  1  a  SECURITY  CUSSFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UL 


Standard  Form  238  -Rev  2*4*9 
Praacrtaad  by  ANSI  Std  9 
299*102 


N$N  7540-01  -290*9900 


TABLE  OF  CONTENTS 


1.  Introduction  .  2 

2.  Formulation  of  the  problem  .  3 

3.  Dyadic  Green's  function 

for  the  unbounded  problem  .  6 


4.  Discussion . 14 

5.  Conclusion .  17 

References  .  19 

Appendix .  21 


ILLUSTRATIONS 


Figure  1.  Geometry  of  the  problem  .  23 

Figure  2.  Rotational  transformations  between  the 
principal  coordinate  system  of  the 


medium  and  that  of  the  problem 


24 


Accesion  For 

NTIS  CRA&J 
DTIC  TAB 
Unannounced  □ 
Justification 


By _ 

Dist;  ibution  / 


Availability  Codes 


Dist 

Avail  and/or 

Special 

A-! 

1 


1.  INTRODUCTION 


The  need  to  study  and  understand  electromagnetic  wave  propagation 
and  scattering  in  anisotropic  media  arises  in  such  diverse  areas 
as  optics,  plasma  physics,  geophysics,  antennas  and  remote  sensing.  In 
some  cases  anisotropy  occurs  as  an  undesirable  parasite;  in  some 
others  it  is  intentionally  introduced  to  serve  some  special  purpose. 
Naturally  occuring  anisotropic  materials  which  are  of  practical  use 
are  very  few  (e.g.,  crystals).  But  thanks  to  the  recent  advances  in 
materials  technology,  a  wide  range  of  anisotropic  materials  with 
diverse  characteristics  are  readily  available  today.  Alongside, 
several  novel  applications  of  anisotropic  materials  are  being 
continuously  invented  [Kobayashi  and  Terakado,  1980;  Paul  and 
Shevgaonkar ,  1981 ] 1 ’ 2 . 

One  of  the  well-established  tools  for  the  analysis  of  electromagnetic 
radiation  and  scattering  problems  is  the  method  of  Green's  functions  [Tai, 

3 

1971]  .  There  have  been  a  few  investigations  in  the  literature  dealing 

with  the  dyadic  Green's  function  for  anisotropic  media.  But  they  are 

either  restricted  to  the  uniaxial  case  [Lee  and  Kong,  1983] 4  or  only 

6  6 

suitable  to  certain  special  applications  [Krowne,  1984a  and  1984b]  ’  . 
Krowne  [1984b]6  has  found  his  DGF  to  be  very  useful  in  problems  such  as 
electromagnetic  field  analyses  of  integrated  circuits  on  anisotropic 
substrates.  Although  in  the  literature  there  are  papers  devoted  to 
electromagnetic  scattering  from  anisotropic  media  of  very  general  type 
[Graglia  and  Uslenghi,  1984]  we  are  particularly  interested  in  analyses 
which  can  offer  some  physical  insight  into  the  scattering  mechanisms  in 
anisotropic  media.  With  these  in  mind  we  seek  to  obtain  the  DGF  of  a  two- 


2 


layer  biaxially  anisotropic  medium  whose  principal  axes  can  have  any 
arbitrary  orientations.  There  are  several  steps  involved  in  our  procedure. 
First  we  obtain  the  DGF  for  the  unbounded  problem.  We  next  formulate  the 
DGF  for  the  two -layer  problem,  express  the  coefficients  in  terms  of  half- 
space  Fresnel  coefficients  and  finally  evaluate  the  various  half- space 
Fresnel  coefficients.  In  this  report  we  confine  ourselves  to  the  task  of 
deriving  the  DGF  for  the  unbounded  case.  Extension  of  this  work  to  the  two- 
layer  problem  will  be  treated  in  another  report. 

The  contents  of  this  report  are  organized  as  follows .  Section  2 
describes  the  formulation  of  the  problem.  In  the  next  section  we 
derive  the  DGF  of  the  unbounded  biaxial  medium.  In  Section  4  we  have 
a  discussion  of  our  result  and  the  report  concludes  with  a  brief 
summary  in  Section  5. 


2.  FORMULATION  OF  THE  PROBLEM 


The  geometry  of  the  two-layer  problem  Is  shown  in  Figure  1.  It 
consists  of  three  regions  -  Region  0  (  z  >  0  )  is  an  isotropic  medium 
with  permittivity  £q,  Region  l(0>z>-d)is  the  anisotropic  medium 
with  permittivity  c  and  Region  2  (  z  <  -d  )  is  an  isotropic  medium 
with  permittivity  c^-  All  three  regions  have  the  same  permeability  /* . 

The  permittivity  of  the  medium  in  Region  1  is  a  second- rank  tensor 
whose  matrix  representation  in  its  principal  coordinates,  take 

the  following  form. 


(1) 
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In  order  to  obtain  the  representation  of  in  the  coordinate  system 

of  our  problem  (x,  y,  z)  we  perform  two  rotational  transformations. 

Let  x" ,  y"  and  z"  be  the  principal  axes  of  the  permittivity  tensor. 
Note  that  the  normal  to  the  boundary  is  not  along  any  of  the  principal 
axes  in  general.  The  first  transformation  (see  Figure  2(a))  is  an 
anticlockwise  rotation  through  an  angle  about  the  x"  axis.  The 
second  transformation  (see  Figure  2(b))  is  an  anticlockwise  rotation 
through  an  angle  ^  about  the  z'  axis.  Thus  the  representation  of  the 
permittivity  tensor  in  the  (x,  y,  z)  coordinate  system  is  given  by  the 
following  symmetric  matrix. 

C11  c12 
c  “  c12  c22 
c13  C23 


13 


'23 


'33 


(2) 


where 


2  2  2  2 
£11  -  cx  cos  i/>2  +  (c  cos  ^  +  ez  sin  sin  if>2 


2  2 

£.,„  -  (-  £  +  €  cos  0,  +  £  sin  i/>. )  sin  cos  i/>0 

12  x  y  1  z  1  2  2 


£13  -  (c2  -  £  )  sin  i>1  cos  ^  sin  ^ 


2  2  2  2 
c22  “  cx  sin  ^2  +  ^ey  cos  ^1  +  cz  Sin  ^  COS  ^2 


£23  -  (£z  ~  £  )  sin  cos  ip1  cos  r/>2 


2  2 

c33  “  cy  Sin  *1  +  £z  COS  *1 


(3a) 


(3b) 


(3c) 


(3d) 


(3e) 


(3f) 


4 


We  have  an  electric  current  source  in  Region  0  away  from  the 
boundary.  We  are  interested  in  finding  the  DGF's  :  Gg^r.r'), 

G^Q(r,r')  and  G2Q(r,r').  The  second  subscript  denotes  the  region 
where  the  source  is  located  and  the  first  subscript  denotes  the  region 
containing  the  observation  point.  In  order  to  obtain  the  above  DGF's 
we  need  to  solve  the  following  equations  [Tai,  1971] 3 : 

V  x  V  x  G00<r.f)  -  J-f.  c0  G00(r,f’>  -li(f-f')  (4) 

V  x  V  x  Gl0(r,r-)  -  1  210(t,t')  -  0  (5) 

V  X  V  x  G20(r,f)  -  Jp  e2  G20<f,r’)  -  0  (6) 

where  w  is  the  angular  frequency.  The  boundary  conditions  associated 
with  these  DGF's  are  given  as  follows: 


z  x  G00<r.r')  -  z  x  G10(r,r')  at  z  -  0  (7a) 

$  x  V  x  G00(r,r')  -  $  x  V  x  G10(r,r')  at  z  -  0  (7b) 

z  x  G10(r,r')  -  z  x  G20(r,r')  at  z  -  -d  (8a) 

z  x  V  x  G1Q(r,r')  -  z  x  V  x  G2Q(r,r')  at  z  -  -d  (8b) 


We  find  this  boundary  value  problem  too  complicated  to  attack  head 
on.  Since  most  of  the  complications  are  due  to  the  anisotropic  medium 
we  first  obtain  the  DGF  of  the  unbounded  anisotropic  medium.  This  is 
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the  problem  we  set  ourselves  to  solve  In  this  report. 


3.  DYADIC  GREEN'S  FUNCTION  FOR  THE  UNBOUNDED  PROBLEM 


The  DGF  of  the  unbounded  problem,  G(r,r'),  satisfies  the  following 
equation: 

V  x  7  x  G(r.r')  -  w2p  7  G(r.f')  -  T  5(f  -  r ' )  (9) 

Although  this  equation  is  very  similar  to  (5),  the  major 
simplification  lies  in  the  absence  of  the  complicated  boundary 
conditions  (7)  and  (8).  The  only  condition  to  be  satisfied  here  is  the 
radiation  condition. 

Instead  of  seeking  solution  to  (9)  directly,  we  find  it  convenient 
to  first  introduce  an  'auxiliary'  dyadic  Green's  function  (ADGF)  3(r,r') 
which  is  linearly  related  to  DGF  (see  (36)).  As  we  shall  see  the 
dyadic  decomposition  is  easier  to  achieve  for  ADGF  which  satisfies  the 
following  equation: 


V  x  V  x  x  3(r,r')  -  w2p  3(r,r')  -  7  6(r  -  r') 


where 


Now  we  define  the  following  Fourier  transform  pair: 


3(k,r')  -  j  d3r  3(r,r')  exp  (-ik*r) 

»(f,r')  -  f  d3k  5(k,r')  exp  (ik»r) 


(10) 

(ID 

(12a) 

(12b) 


Then  in  the  Fourier  domain,  (10)  reduces  to  the  following  algebraic 
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equation. 


k  x  k  x  x  3(k,r')  +  u>  ft  d(k,r')  -  -  c  exp  (-ik*r’) 


The  solution  to  (13)  is  readily  obtained  as 


(13) 


3(k,r')  -  -[k2;  +  4i]  1  “ 


c  exp  (-i  k*r' ) 


(14) 


where  the  matrix  representation  of  k  in  the  (x,  y,  z)  coordinate 


system  is  given  as 

k  - 


0  -k  k 
z  y 


k  0 

z 

-k  k 

y  * 


-k 


(15) 


Substituting  (14)  in  (12b) ,  we  get 


3(t,t') 


“  (2*)s  I 


3—  J  adJ  (  k2  x  +  A  l) 


d^k 


^  det  [  k^  x  +  l] 


e  exp 


[ik.(r-r')]  | 


(16) 


r  -2  -  2  -  1 

Factorizing  det  Ik  ^  +  u  /i  I  I  ,  we  obtain 

det  [  k2  x  +  A  T  ]  -  A  det  x  (  kz’k*U  ) 


(  Vkf  )  {  k,-k^  )  (  k.-k“  )  (17) 


,bu 

z  z 


bd 
z  z 


Here 


au 

z 

ad 

z 


-  |  (-  A/2  ±  *  ±  SD  ) 


(18a) 


7 


f  ,  bu 

k 

z 


k 


bd 


-  \  (-  A/2  ±  *  ±  «  ) 


(18b) 


where  A,  S,  S  and  J£  are  given  in  the  Appendix.  The  superscripts  a  and 

b  refer  to  the  two  characteristic  waves  that  exist  in  a  biaxially 

anisotropic  medium.  Hereafter  we  shall  denote  them  as  a-wave  and  fa- 

wave.  The  superscripts  u  and  d  denote  the  upward  and  downward 

propagating  waves,  respectively.  Note  that  Re  kaU  and  Re  k^u  are 

z  z 

positive  while  Re  kaC*  and  Re  k^  are  negative.  These  are  in  agreement 

with  our  convention  for  upward  and  downward  travelling  waves. 

Furthermore,  we  assume  that  the  medium  is  slightly  lossy  so  that 

Re  kfi  ►  1m  kfi,  C  -  a  or  b  and  i  -  u  or  d.  Hereafter  in  this  paper  we 
z  z 

follow  the  above  notation  for  f  and  i. 

From  (16)  and  (17)  we  obtain 


--  (2^ 


ik  •(p-p')  ik  (z-z '  ) 
P  Z 

e  r  e 


adj  (  k* 


X  +  o>  p 


i) 


(kz-kzU^kz-kzd)(kz-kzU)(kz'kzd> 


(19) 


We  now  perform  the  integration  in  (19)  over  k^.  The  integration  is 

from  -<*>  to  -H»  along  the  real  axis.  There  are  four  poles  in  the  complex 

k  -  plane,  viz.,  at  kau,  kaC*,  k^*U  and  k*5^.  The  poles  at  kau  and  k 
z  z  z  z  z  z  z 

lie  in  the  first  quadrant  while  the  poles  at  ka<*  and  k^  lie  in  the 

z  z 

third  quadrant.  For  z  >  z'  we  deform  the  integration  contour  in  the 

au  .  ,  bu 

upper  half-plane  and  this  contour  encloses  the  poles  at  k  and  k  . 
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Using  Cauchy's  theorem  the  following  result  is  readily  obtained. 
For  z  >  z' 


3(r,r')  -  - 


c  c  c 

i  xyz 

(2jt) 2  a7?  e,. 


where 


r  , au  au  A+  A+ 

J  _ X  g  a  a _ 

^  <kauVd)(kau.kbu)(kau-kbd) 

z  z  z  z  z  z 


exp  [ika»(r-r')j 


,bu  bu  A+  A+ 

_ A  g  4  b _ 

<kb“.ka“)(kbu.kad)(kbu.kbd) 

z  z  z  z  z  z 


exp  [ikb»(r-r')j  j 


A+ 

a  — 


£a  a 
k  x  o- 


|kax 


.  a  f\ 

k  x  o2 

,  ^a  a  . 
|k  x  o2| 


(21a) 


^b  .^b  A  .  ^b  .^b  A  . 

k  x(k  x  o^)  k  x(k  x  o 2> 

,  *b  A  .  +  .  ^b  A  . 

|k  x  Ql|  |k  x  o2| 


(21b) 


£  k  +  $  kru 
P  z 


f  -  a  or  b 


(kr)2 

«2/i  -  (kr)2 


gf"  -  (k^4  ( xf2u  i2  +  [»•<.-  (kci2  r 
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o>  •  o> 


f  (  kf  x  ox  ) • (  kr  x  o2  )  ^  1/2 

^  |kf  x  o^  |kf  x  o2|  f 


kf  -tf/k* 


k^  •»  |  k^ 


61  " 


£z<V  £x>  V2 

c  (c  -  c  ) 
y  z  x  J 


cx(£z~  V  ^ 

c  (c  -  c  ) 
y  z  x  J 


fu  ^  2  ^fu  J  2  ,Cu  „  ,  ,fu  .{ 

X1  “  X11  sin  ^  +  x22  COS  ^  "  2  xi2  sin  ^  cos  * 


(  Xn  -  X22  )  cos  sin  ^ U  cos  ^ 


,  -fu  r  2  .fu  2  .fu  ■» 

+  x^2  cos  6  [  sin  p  -  cos  <f>  J 


+  (  X9o  cos  -  x,,  sin  <f>*U  )  sin 


(26a) 


(26b) 


"  X  (  ±  cos  \f>2  +  g2  sin  ^  sin  +2  ) 

+  y  (  *  gj^  sin  +  g2  sin  ^  cos  V>2  )  +  *  cos  ^  (27) 


(28a) 


(28b) 


(29a) 


(29b) 
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Xjj  is  the  ij  element  of  the  matrix  x  I  the  angles  ($^U,  ] 

Af 

denote  the  direction  of  k* .  For  definiteness  we  have  assumed  in  the  above 

equations  that  c  <  c  <  c 
x  y  z 

For  z  <  z'  we  deform  the  integration  contour  of  (19)  in  the  lower 

half-plane  thus  enclosing  the  poles  at  kad  and  kbd.  Once  again,  use  of 

z  z 

Cauchy's  theorem  leads  to  the  following  result. 


For  z  <  z' 


3<r,r')  - 


A _ JLXJL  (Vi; 

?*>2  c 33  J  d  % 


r  .ad  ad  a-  a- 

J  _ X  g  a  a _ 

1  „ad  .  au.  ..  ad  ,  bu.  ..  ad  ,  bd. 

1  (kz  -kz  )<kz  -kz  Hkz  -k2  ) 


exp  [ixa*(r-r')l 


xbd  gbd  r  r 

<kbd-k‘uMkbd-kad>(kbd.kbu> 
z  z  z  z  z  z 


[i*b* (r-r' )]  | 


where 


,xa  A 
k  x  o. 


,  Aa  A 

I*  x  0l( 


Aa  a 
K  X  02 

.Aa  a  . 

I*  X  o2l 


Ola) 


Ab  Ab  A  . 
K  X(«  X  O^) 

i  Ab  A  . 

\k  x  0l| 


Ab  Ab  A  . 
K  X(K  X  O2) 

iAb  a  , 
I*  x  o2| 


(31b) 
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— f  A  ,  .  A  .  f  d  _  , 

#e  -  p  k  +  z  k*  ,  c-aorb 

p  z 

(32a) 

1 

(32b) 

Afd  -  xfu  {  u  -  d  } 

(33) 

gfd  -  gfU  {  u  -  d  } 

(34) 

hrd  -  hru  {  u  -  d  } 

(35) 

Now  that  we  have  derived  the  ADGF  3(r,r')t  it  is  straightforward  to 
obtain  the  DGF  G(r,r').  On  post -multiplying  (10)  by  x  and  comparing  it 
with  (9),  we  deduce  that  the  DGF  is  linearly  related  to  the  ADGF  as 


G(r.r')  -  x  3(r,r')  x  (36) 

Thus  the  DGF  for  the  unbounded  problem  is  finally  written  as  follows . 
For  z  >  z' 


G(r , r ' ) 


i 

(2w)  2 


c  c  e 
x  y  z 


w2p  e 


33 


I 

'  1 


au  (au  au  A+  A+ 
v  A  g  a  a 

<kau-ka<i)  (kau-kbu)  (kau-kbd) 
z  z  z  z  z  z 


exp  ^ika*(r-r')j 


+ 


^/>u  ^bu  ^bu  £+  £+ 


(kbu-kau)(kbu-kad)(kbu-kbd) 
z  z  z  z  z  z 


and  for  z  <  z' 


(37) 
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{ 


ad  ad  ad  A-  A- 
v  A  g  a  a _ 


exp  [i*a*(r-r')j 


4d  .bd  bd  a-  a- 
v  A  g  d  b  r  -b _ l  1 

+  ,.bd  .«uw.bd  .adwlbd  ,bu  exP  liK  f 

(k  -k  )  (k  -k  )<k  -k  )  > 


where 


A+ 

a 


A- 

a 


fi+ 

r 


(  uau)'1/2 


( >.**)' 1/2 


*  3-/2 

*  J 


( *w)‘l/2 


X 


X 


X 


X 


A+ 

a 


A- 


a 


&+ 


r 


tu 

1/ 


l  -  a  or  4 


Id 

v 


A —  _o  A_ 

«■  .  x‘  •  *■ 


t  -  a  or  4 


(38) 

(39a) 

(39b) 

(40a) 

(40b) 

(41a) 

(41b) 


13 


4.  DISCUSSION 


Having  now  derived  the  DGF  of  the  unbounded  biaxially  anisotropic 
medium,  it  will  be  instructive  to  study  some  of  its  characteristics. 
First  we  compare  our  DGF  with  those  of  a  corresponding  uniaxial 
problem  [Lee  and  Kong,  1983]4.  We  notice  that  our  result  is  identical 
in  structure  to  theirs.  The  most  important  difference  is  that  in  the 
place  of  ordinary  and  extraordinary  waves  we  have  a-  and  b-  waves. 

In  order  to  study  the  characteristics  of  a-  and  b-  waves  we  first 
examine  the  dispersion  equation 


det  [  x  + 


u2 n  T  ] 


When  we  solve  (42)  for  k  we  have  found  in  (18)  that  there  are  four 

z 

2 

solutions.  However,  when  we  solve  (42)  for  k  we  find  that  it  has  two 

a  2  b  2 

solutions,  viz.,  (k  )  and  (k  )  ,  where 


a**)2 


r  r  2  -  n  1/2  -| 

—  t  +  t  -  4  d 
d  L  a  v  a  aJ  J 


(43a) 


<*V 


-  Ac  -  Af 
t^  ■  tr  x  ”  k  •  X  •  k 


\  -( 


c  c  c 
x  y  z 


Ar  -  Af 

ks*  £  .  k^ 


(43b) 


(44a) 


(44b) 


Si  b 

It  is  clear  from  (43)  that  the  propagation  constants  k  and  k 
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functions  of  the  direction  of  propagation.  This  implies  that  both  a- 

and  b-  waves  belong  to  the  extraordinary  category. 

But  when  we  let  c  -  c  and  -  0* ,  (A3)  becomes 
x  y  i. 

(ka)2  -  ex  (45a) 

<kV  -  o>2m  rz  +  (k*U)2  (  1  -  rzAx  )  (A5b) 

We  now  see  that  k^  is  dependent  on  the  direction  of  propagation  while 
ka  is  not.  Thus  in  this  special  situation  the  a-wave  is  an  ordinary 
wave  and  the  b-wave  is  an  extraordinary  wave.  This  is  a  well-known 
result  of  a  uniaxial  medium. 

Further  we  note  from  (18)  that  the  magnitudes  of  the  z- component  of 

the  propagation  vectors  are  different  for  upward  and  downward 

travelling  waves,  i.e.,  |k^Uj  +  |k^d) .  This  is  true  in  a  general  case. 

z  z 

But  there  are  situations  where  these  magnitudes  become  equal.  One  such 
situation  occurs  when  the  principal  axes  of  the  medium  coincide  with 
(x,  y,  z),  i.e.,  when  ^  ^  •  In  we  note  that  A  -  R  - 

0.  Thus  it  is  clear  from  (18)  that  z- components  of  the  propagation 
vectors  of  the  upward  and  downv^iu  travelling  waves  are  equal  in 
magnitude  but  opposite  in  sign.  This  result  is  intuitively  satisfying. 

We  next  turn  our  attention  to  the  polarization  vectors.  We  note 
that  our  a  and  %  correspond  to  and  in  Njoku  [1983] 8  and  Kong 
[1975]°.  They  have  denoted  the  characteristic  waves  as  type  I  and  type 
II  waves.  We  further  notice  from  (21)  that  for  a  specified  propagation 

Aa  Ab  A*t* 

direction  (k  -  k  ) ,  a  •  b  -0.  A  similar  result  has  been  reached  by 
Kong  [ 1975 ] °  and  Njoku  [1983]*  by  a  different  method.  Also  from  (21)  and 
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(40a)  we  obtain  that  a+«  6+-  0. 

Although  a+  and  &+  are  orthogonal,  it  is  interesting  to  note  that 
a+  and  £+are  not,  in  general,  orthogonal  to  each  other.  For  we  find 
from  (39)  and  (40)  that  a+»  6+-  0  only  if  k  x  o^  -  ±  (k  x  o^).  This 

implies  that  a+  and  are  orthogonal  only  in  the  special  situation 

A  A  A 

where  k  either  lies  in  the  plane  bisecting  the  angle  between  o^  and  o ^ 

or  is  normal  to  the  bisecting  plane. 

On  the  other  hand,  in  the  case  of  a  uniaxial  medium  the 

corresponding  vectors  are  always  orthogonal  [Kong,  1975;  Lee  and  Kong, 

1983]  ’  .  In  order  to  clarify  this  apparent  difference  we  proceed  as 

follows.  First  we  take  the  uniaxial  limit  of  our  results.  In  this 

limit  we  note  from  (45)  that  (ka)^  -  (k°)^  and  (k'3)^  -  (ke)^.  The 

superscripts  o  an  e  refer  to  the  ordinary  and  extraordinary  waves  as 

in  Lee  and  Kong  [1983] 4 .  Further  (27)  and  (21)  reduces  to  the  following 

limiting  values. 

o^  -  c>2  -  y  sin  ^  +  z  cos  ^  (46) 


A+ 

a  — 


^o  A 
k  X  o^ 

,  f'o  A  , 
|k  x  OjJ 


(47a) 


^e  ,(^e  A  . 

k  x(k  xo^) 

.  ^e  ,^e  A  .  . 

|k  x(k  xo^) | 


(47b) 


Also  in  this  uniaxial  limit 


A+  A+ 
a  -  a 


Since  we  have  seen  earlier  that  a+»  6+-  0,  it  follows  that  a+*  6+-  0 
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always  In  the  uniaxial  limit.  Thus  all  the  properties  of  the  uniaxial 
medium  fit  logically  well  within  the  framework  of  the  biaxial  medium. 

This  completes  the  first  part  in  our  derivation  of  the  DGF's  of 
the  two -layer  biaxially  anisotropic  medium.  In  the  second  part  [see  also 
Mudaliar  and  Lee,  1991] 10  we  extend  this  work  to  the  corresponding  two- 
layer  problem.  The  most  important  merit  of  this  DGF  lies  in  the 
repesentation.  Similar  representation  [Lee  and  Kong,  1983] *  has  served 
us  well  in  studying  and  analyzing  various  scattering  problems  in 
uniaxially  anisotropic  random  media  [Lee  and  Kong,  1985;  Mudaliar  and 
Lee,  1990] 11 ’ 12 .  It  is  expected  that  the  DGF  derived  here  will  serve 
similar  purposes  in  radiation  and  scattering  problems  involving 
biaxially  anisotropic  media. 

5.  CONCLUSION 


Due  to  numerous  applications  in  various  areas  anisotropic  medium  is 
demanding  increasing  attention  in  recent  years.  Since  dyadic  Green's 
function  (DGF)  is  a  basic  tool  for  analysis,  we  seek  to  derive  the  DGF 
for  a  biaxially  anisotropic  two -layer  medium.  To  this  end  we  derived  in 
this  paper  the  DGF  for  an  unbounded  anisotropic  medium  whose  principal 
axes  are  arbitrarily  oriented.  The  Fourier  transform  method  was  used  to 
derive  an  auxiliary  dyadic  Green's  function  (ADGF) .  The  DGF  was  then 
obtained  by  performing  a  linear  transformation  on  the  ADGF.  Some 
properties  of  the  characteristic  waves,  viz.,  the  a-waves  and  the  b- 
waves,  were  studied.  Our  representation  of  the  DGF  demands  particular 
attention.  Besides  its  compact  structure  it  also  has  physical 
correspondence  to  the  two  characteristic  waves  in  the  medium.  In  the 
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uniaxial  limit  our  results  agree  well  with  those  in  the  literature. 
Extension  of  this  work  to  the  two-layer  problem  will  be  discussed  in 
another  report. 
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APPENDIX 


A  ■  2<kx  *13  +  ky  c23)/c33 
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